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SUMMARY 

i 

A typical low-thrust Mars orhiter mission using the Snap-8 power- 
generating system has been studied from a guidance viewpoint. The 
mission trajectory was divided into the three characteristic phases, 
namely, the outward escape spiral, the heliocentric transfer, and the 
inward capture spiral. The sensitivity of the final trajectory state 
(velocity, position, and mass) due to errors in the initial trajectory 
state and to thrust-vector errors was determined for each phase of the 
mission. The analysis and numerical solutions are based on the methods 
of linear perturbations and adjoint functions. 

During the escape -spiral phase, any of the following minimum errors, 
if left uncorrected, is sufficient to cause the escape direction to be 
180° out of phase: a 0.1-percent error in either initial orbital alti- 

tude or initial vehicle mass, a 0.1-percent error in the thrust magnitude 
acting over the entire trajectory, or a 2.8° error in the thrust angle 
acting over the entire trajectory. 

During the heliocentric transfer, final position and velocity errors 
of the order of several hundred thousand kilometers and 10 meters per 
second will result r from either an escape -direction error of l/2°, a 
thrust -magnitude error of 1 percent, or a thrust-angle error of l/2°. 

During the capture-spiral phase, if the initial velocity vector were 
in error by only 1 meter per second and 1 milliradian, and if the nominal 
thrust program were followed exactly, the vehicle would spiral down into 
the Martian surface. 

In general, low -thrust trajectories are highly sensitive to errors. 
Open-loop trajectory control is thus completely out of the question; that 
is, repetitive trajectory determination and corrective guidance maneuvers 
are required to ensure a successful mission. 
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INTRODUCTION 

In the past several years, theoretical and experimental analyses 
have shown that electric rocket systems are practical and offer high 
performance in advanced space missions. First-generation spacecraft 
flights are scheduled to begin around 1965 with an ion engine powered 
by the Snap-8 nuclear turboelectric system currently under development. 

A feasible mission that has been proposed for the 60-kilowatt version 
of Snap-8 is a scientifically instrumented Mars orbiter. Such a mission 
would begin in a low-altitude Earth-satellite orbit and end in a similar 
orbit about Mars-. Since the vehicle would penetrate a large portion of 
the near-Martian space as it slowly transverses the capture -spiral tra- 
jectory, an accurate survey of radiation and surface conditions may be 
obtained and transmitted back to Earth. 

The question of how accurately the vehicle can be guided to its 
given target is of great consequence to the success of the mission. As 
yet, little attention has been given to the problem of low-thrust guid- 
ance in the presence of random or systematic perturbations arising from 
such error sources as thrust-vector control, navigational measurements, 
and an approximate system model. In a preliminary analysis of this 
subject (ref. l), linear perturbation theory and the method of adjoint 
functions were used to derive the fundamental guidance equation for 
low-thrust trajectories. This equation provides a means of studying 
the effect of error perturbations and corrective guidance perturbations. 
The specific problem treated in reference 1 is an error analysis of the 
heliocentric-transfer trajectory between Earth and Mars covering a wide 
range of propulsion characteristics and transfer times. This report is 
an investigation of the guidance problem for a typical Mars orbiter 
mission considering, in turn, the escape-spiral phase, the heliocentric- 
transfer phase, and the capture-spiral phase. The purpose of this in- 
vestigation is to determine the sensitivity of the final trajectory state 
for each phase to initial trajectory state errors and thrust vector 
errors. It should be emphasized, that this study is a fixed-time anal- 
ysis; that is, perturbations refer to the difference between the actual 
and reference trajectories at a particular instant of time. The choice 
of time as the independent variable is computationally convenient and 
is a good criterion of comparison for the problem at hand. 


SYMBOLS 

a semimajor axis, m 

c jet velocity, m/sec 

e eccentricity 
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F thrust force, newtons 

function of thrust vector perturbations (eq_. (19)) 

GM gravitational constant of central attracting body, m 3 /sec 2 

m vehicle mass, kg 

p semilatus rectum, m 

r radial position, m 

r a apogee distance, m 

r^ perigee distance, m 

t time, sec or days 

u radial velocity, m/sec 

V total velocity, m/ sec 

W matrix of weighting functions 

w thrust-error weighting function 

x general trajectory state variable 

y general thrust -vector perturbation variable 

0 thrust angle, radians or deg 

A( ) large variation from reference quantity ( ) 

&( ) small variation from reference quantity ( ) 

A matrix of sensitivity coefficients 

A error sensitivity coefficient, variable in adjoint equations (8) 

to (12) 

cp angular position, radians 

cd angular velocity, radians/sec 

• scalar product of vectors 
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Subscripts: 

c beginning of coast period 

des design thrust level 

E Earth 

f final value 

i j. general indexes 

M Mars 

n values at t n - beginning of control interval 

0 initial value 
Superscripts: 

T transpose of matrix 

1 second-order perturbation 
time derivative 

vector 


ANALYSIS 

Vehicle and Mission Characteristics 

In this report an investigation of the trajectory sensitivity and 
control problem for a particular Earth-Mars instrumented space probe 
utilizing the Snap-8 electric propulsion system is presented. Rather 
than proceeding directly to the analysis, it is appropriate here to 
present the propulsion characteristics of the vehicle that have been 
assumed and to discuss briefly the mission and trajectory characteristics. 

Vehicle-propulsion characteristics (unpublished NASA data) . - The 
Snap-8 vehicle will be placed into a nearly circular, low-altitude 
satellite orbit utilizing a multistage! chemical booster. Typically, 
an initial gross weight of 4080 kilograms (9000 lb) is placed in a 
927-kilometer (500-naut.-mile) orbit. The design electric power and 
specific powerplant mass are 60 kilowatts and 22.7 kilograms per kilo- 
watt, respectively, and the resulting powerplant - initial-gross-weight 
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ratio is l/3. With a propellant utilization of 90 percent and a power 
efficiency of 76.2 percent assumed, a design specific impulse of 3600 
seconds is obtained. The available thrust force is 2.32 newtons 
(0.524 lb), and the propellant mass rate is 5.68 kilograms per day. 

Mission characteristics . - The mission under consideration is a 
415.4-day, one-way trip to Mars beginning in a 927 -kilometer Earth- 
satellite orbit and terminating in a similar satellite orbit about Mars. 
Figure 1 is a schematic illustration of the overall interplanetary- 
transfer maneuver. It is convenient to discuss the mission in terms of 
three distinct phases, namely, the escape spiral, the heliocentric trans- 
fer, and the capture spiral. 

(1) During the first 139 days of the trip the vehicle increases its 
energy relative to Earth according to a continuous -power tangential - 
thrust program. Solar gravitational effects are neglected during this 
phase. 

(2) At the hyperbolic-escape condition, the gravitational effect of 
the Earth is neglected, and the vehicle initiates a 223-day heliocentric 
transfer employing an optimum constant -thrust program. The power-coast- 
power sequence is 27.7, 160, 35.3 days, respectively. The heliocentric- 
transfer angle is 159°. 

(3) At 362 days the vehicle is in the vicinity of Mars and initiates 
the energy-decreasing capture spiral, again employing a continuous-power- 
tangential -thrust program. After 53.4 days, the 927-kilometer satellite 
orbit has been established. Solar gravitational effects are neglected 
during this phase. 

The overall trajectory has been pieced together from a series of 
two-body solutions in order to simplify the analysis. This procedure 
has been commonly used by investigators of the low-thrust interplanetary 
mission. In most cases, the approach taken is to terminate the escape 
spiral when parabolic escape energy has been attained. After this, the 
vehicle is assumed to be moving in the Earth's orbit with Earth's orbital 
velocity relative to the sun, and the Earth's gravitational effect is 
then neglected. A similar procedure is used at Mars. Frequently, the 
parabolic energy condition occurs at a position well within the planet's 
sphere of influence. The essential difference in the approach taken in 
this analysis is an extension of the escape- and capture- spiral phases 
past the parabolic energy condition. Transformation to or from the 
heliocentric system is arbitrarily made at a point where the solar and 
the planetary gravitational effects are approximately equal. The main 
reason for extension of the spiral phases is that in the vicinity of 
the hyperbolic escape or capture condition the motion of the vehicle is 
approximately along the asymptote of a hyperbola, which is a useful 
target criterion for guidance purposes. 
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Figure 2 describes the escape-spiral trajectory in detail, where 
the last of 500 revolutions about Earth is plotted. Parabolic escape 
energy is reached at a distance of about 100 Earth radii and at a time 
of 125.5 days after launch. At a distance of 298 Earth radii, the escape 
phase is terminated, and, at this point, the relative velocity of the 
vehicle is 1577 meters per second. The direction of the velocity vector 
was arbitrarily chosen to be parallel to Earth's velocity. The perpen- 
dicular distance from the center of the Earth to the hyperbolic asymptote 
is referred to as the asymptotic displacement and, in this case, is 
0.675X109 meters.. 

As previously mentioned, the heliocentric-transfer trajectory is 
of the optimum constant -thrust type; that is, the thrust program employed 
was one that minimized the propellant expenditure subject to the follow- 
ing constraints: (l) The thrust must be equal to its maximum design 

value or equal to zero, and (2) the transfer is to take place in a spec- 
ified time, namely, 223 days. Details of the optimization procedure 
are fully presented in reference 2. 

A plot of the capture -spiral trajectory is shown in figure 3. Ini- 
tially, the vehicle ' s relative position and velocity are 418 Mars radii 
and 1454 meters per second, respectively. The direction of relative 
velocity is parallel but opposite to Mars' orbital velocity. Thus, the 
planet is catching up with the vehicle. The asymptote of the capture 
hyperbola is displaced by 0.209X1C>9 meters. With the use of reverse 
tangential thrust, a 927-kilometer satellite orbit is achieved during 
the 109^ spiral turn. 


Trajectory Sensitivity Analysis 

The motion of the vehicle during each phase of the mission is de- 
scribed by a set of nonlinear differential equations. Consider that a 
solution of such a set corresponding to a specified thrust program and 
satisfying prescribed boundary conditions has been obtained and is 
termed the reference solution. If the major assumption is made that 
perturbations (e.g., guidance errors) are sufficiently small so that 
the actual vehicle trajectory does not vary significantly from the 
reference trajectory, it is possible to study these variations and the 
required corrective maneuvers by linear perturbation techniques. The 
perturbed differential equations of motion represent a linear system with 
time-varying coefficients, which are determined from the known reference 
trajectory. The general solution of the perturbed equations is best 
facilitated by the method of adjoint functions as suggested in refer- 
ence 3. Application of the adjoint method to various trajectory sensi- 
tivity and control problems is found extensively throughout the litera- 
ture, for example, references 1, 4, and 5. 


7 


Equations of motion . - During each phase of the mission, the vehicle 
is ass ume d to travel in a vacuum under the influence of an inverse -square 
central gravitational field in addition to its own thrust acceleration. 

A two-dimensional geometry is used, and the vehicle motion is described 
in a rotating polar coordinate system centered at the appropriate central 
body. With reference to figure 4, the differential equations of motion 
to be satisfied along the flight path are given as 


>2 _ GM 

1? 

+ — sin (3 

(1) 

•j.2 

m 


-2uo> j 

F 

— cos p 

(2) 

r 

rm 


r = 

u 

(3) 

9 = 

U) 

(4) 

m = • 

. 1 

(5) 


Differentiation with respect to the independent variable time is denoted 
by a superscribed dot. The state variables of interest are radial ve- 
locity u, angular velocity u>, radial position r, angular position 9, 
and mass m. The propulsion or control variables are the thrust magnitude 
F and the thrust angle (3, which is measured with respect to the local 
horizontal (circumferential) direction. The effective jet velocity c 
is equal to the product of specific impulse and a conversion factor and 
is considered constant in this analysis. 

The fact that the angular position 9 does not appear in the non- 
linear equations of motion is very significant in that large variations 
in 9 may be admitted without invalidating the linear perturbation 
analysis. The possibility of very large variations in this quantity, 
relative to 2 jt radians, during the escape-spiral phase is evident, since 
the vehicle makes several hundred revolutions around Earth in the process 
of escaping. 

Fundamental guidance equation . - The derivation of the linearized 
equations of motion and their solution by adjoint methods are reported 
in reference 1. The development needed for the present study is presented 
in appendix A. In the following discussion, the symbol & is used to 
represent small variations from reference quantities, that is, by def- 
inition 


Sx(t) = x^t) - x*(t) 
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where x is any state or control variable, and the asterisk denotes 
the reference value. 

The fundamental guidance equation expresses the variation in state 
variables at the final reference time tp in terms of the variation in 
state variables at some time t along the path and in terms of the 
integrated effect of thrust-vector variations during the interval 
(t, tf ) . This equation is written with the use of matrix notation for 
conciseness, as 



where A is a 5 by 5 matrix of sensitivity coefficients and W is a 
5 by 2 matrix of thrust sensitivity or weighting functions. By 
definition 



The elements Apj and w^j are determined from a set of adjoint dif- 
ferential equations discussed in the next section. 

The fundamental guidance equation has two basic interpretations. 
First, the sensitivity of final conditions due to a variety of error 
sources may be determined, and the analyst thereby provided with in- 
formation regarding the navigational and control-accuracy requirements. 
For example, if the final variations have some prescribed tolerance 
level, (&F, 50) may be regarded as constant bias or random variations 
acting over the path, and their tolerance level may be computed. Also, 
if at time t the variations 5u, 5<n, and so forth are measured, the 
uncertainty in these variations due to instrumentation errors causes 
uncertainty in the knowledge of final conditions that may be determined 
from equation (6). Thus, something might be inferred about the measure- 
ment accuracy requirements . The second interpretation of the equation 
has to do with the formulation of an active guidance technique. With 
the assumption that corrective maneuvers are required to null some or 
all final variations, which have been determined by measurements Su(t), 


) 


9 


5cu(t), and so forth, the necessary corrective-thrust program is implic- 
itly contained in the integral terms of equation (6). The form of 
(5F, 8(3) will depend on the particular control criteria and constraints 
imposed . 

In this report the major concern in the evaluation of the effect of 
errors in initial conditions and thrust vector on the final-state vari- 
ables, that is, the evaluation of the right side of equation (6) when 
t = t Q . For example, if the error in radial velocity at tf is of 
interest equation (6) yields 

8u f = (A-jj_8u + A 12 6 co + A 13 8r + A-^Scp + A 15 8m)t 0 


+ 



(wn&F + w 12 &3)cit 


From the form of this equation, the sensitivity coefficients may be 
interpreted as partial derivatives j for example, 


C*Uf 

A H (t o> s ^ 


Also, if &F and 83 are constant along the path, the following partial 
derivatives could be defined as 


ST = 


The extension of these definitions to the remaining variables follows 
directly. 

Adjoint equations . - The sensitivity coefficients Aij(t) are deter- 
mined from the solution of the following adjoint differential equations 
(see appendix A) : 




>ul = \L2 - ^i5 


( 8 ) 
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A. 


12 = " (2rcD)A iX + (2|)A i2 - A i4 


+ rf)»n - (|F - ^>13 


A i4 = 0 


is - il ^ 


rim- 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 


The time-varying coefficients in parenthesis are known from the particular 
reference trajectory of interest. These equations were integrated 
backward in time from tf. The initial conditions specified at tp have 
the form 


Aii(tf) = 1) A i j(t f ) =0 j ^ i 


that is, 


A(tf) = l(unit matrix) 


(13) 


The sensitivity coefficients and p are continuous functions of 
time even if a discontinuity in F is admitted, as it is during nominal 
coasting periods (see ref. 2). From equations (ll) and (12) it is noted 
that A 4 4 is a constant and that Aqs is constant during coasting 
periods (F = 0). When boundary conditions are specified as in equa- 
tion (13), the fourth column and the fifth row of the matrix A are 
time invariant and may be given directly without numerical integration: 


( A 14; a 24; A 34; a 44> a 54) = (0,0, 0,1,0) 

( A 51> A 52; A53, A54, A55) = (0,0, 0,0,1) 

J 

The elements of the matrix W are, from appendix A, 


A il 

Wii = -y— sin p + — cos p 


m 


F 


A i2 

rm 


A, 


15 

c 


A. 


w i2 = m\ Ai l cos P 


12 


sin p 


(14) 


(15) 

(16) 
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On substitution from the second of equations (14), w,-^ and w^g sim- 
plify to 

W 51 = - b x 52 = 0 (17) 

The thrust-angle weighting function w^g is continuous only during the 

nominal powered periods and is zero during the nominal coast period 
(F = 0) . A discontinuity exists at the initiation and the termination 
of the coast period. This result is as it should he, since thrust-angle 
variations have no physical meaning when thrust is shut off. 

Honlinear thrust-vector perturbations . - In the derivation of the 
fundamental guidance equation (eq. ( 6 ) ) , it is assumed herein that thrust- 
vector perturbations are sufficiently small so that only first-order or 
linear variations need be considered. If, however, large variations in 
thrust magnitude or direction are allowed, the linear form of the inte- 
grand in equation (6) must be modified; that is, 


W 


€>F 

5(3 


[f ± (AF, A(3 ) ] 


i = 1, 
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where AF and A3 may be considered large variations. The need for 
such modification arises when consideration is made of the problem of 
guidance maneuvers (corrective thrust programs). For example, if the 
guidance maneuver requires that thrust be cut off for a specified time 
during the nominal powered period or that thrust be turned on during 
a nominal coast period, than AF = +F, which is certainly not small. 
Another reason for deriving the exact fq functions is that a check on 
the conclusions drawn from the linear -sensitivity analysis is available 
when the first-order sensitivities appear to be negligible. 

The perturbation analysis (see appendix A) involves two essentially 
unrelated types of linearization, that is, linearization of the state 
variables and linearization of the control variables. This separation 
is evident in the results of equation (6). The validity of the trajectory 
perturbation solution, however, depends only on the validity of the 
linearization of state variables.' Wote that the adjoint solution A(t) 
is independent of control -variable linearizations. Consequently, as 
long as large thrust-vector perturbations do not result in large trajec- 
tory perturbations, it is necessary to modify only the integrand terms 
in equation (6) . From appendix A the result is 

Sx(tf) = A(t)5x(t) + 



[fildt 


(18) 
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where 


A il 


f.j_ = -jjp [ (F + AF)sin(p + Ap) - F sin p] 


A. 


12 


A i5 


+ — [ (F + AF )'Cos(p + Ap) - F cos p] [ (F + AF) - F] 


After expansion of the trigonometric functions and with the use of 
equations (15) and (16), the alternative form is 


fi = w i:L 


+ 


A i5 


!AF cos Ap + 


'^12 
V F 


AS? sin Ap - 


A. 


15 


AF 


A. 


+ Fiwj_q + — ^)(cos Ap - l) + (wj_2) s i n AP 


(19) 


When AF and AP are small, AF « 5F, sin AP » Sp, and cos AP 1. 
Equation (19) thus degenerates into the approximate linear form derived 
previously : 


w il 5F 


+ w. 2 Sp 


In the previous section it is noted that wq2 I s discontinuous 

at the beginning and at the end of a coast period (F = 0) and is zero 
during this period. The quantity w-^/F, which appears in equation (19), 
however, is always continuous and, in general, is nonzero, as can be 
seen from equation (16 ) . 

Several special cases of equation (19), used later in this report 
are as follows: 


F = 0; AF = F^ es ; Ap arbitrary 


F i F des 


,w il + 


A i5 


cos Ap + 


w. 


12 
F ) 


sin Ap - 


A i5 


Ap = 0; AF arbitrary 


fi = w ±1 AF 


F = Fq es ; AF = SF « F; Ap = sp, small but second-order effects 
significant; sin AP « 6P; cos AP - 1 « - ^ &P^ 


( 20 ) 


( 21 ) 
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F-3 pc , / Aj et\ 

fi = w u&F + w i2 S J3 - -g^hfu + BP 2 ( 22 ) 

RESULTS MD DISCUSSION 

The overall low-thrust trajectory for the 415.4-day Mars orhiter 
mission has been pieced together from a series of two-body solutions. 

The equations needed to calculate the reference trajectories were pro- 
grammed for an IBM 704 digital computer. Numerical integration was 
performed by the Runge-Kutta technique with an automatic step-size con- 
trol to limit truncation error. The computer program was developed for 
the study presented in reference 2, and only minor modifications were 
required to adapt it to this investigation. Basically, this modification 
involved increasing the integration loop by the addition of the adjoint 
equations (8) to (12). 

In the following sections, each phase of the mission is considered 
separately, and the results show the sensitivity of the final trajectory 
state due to errors in the initial trajectory state and thrust -magnitude 
and -direction errors. Two examples of simple corrective guidance are 
presented. 


Trajectory Sensitivity - Escape -Spiral Phase 

Reference trajectory . - Characteristic parameters of the escape- 
spiral trajectory are given in figure 5. The thrust angle (3, which 
corresponds to the tangential-thrust program, is plotted in figure 5(a). 
Thrust direction is seen to be within 2° of the horizontal throughout 
the first 100 days, and thus a circumferential -thrust program would 
result in approximately the same trajectory up to this point as the 
tangential program. Aside from the fact that tangential thrust is very 
efficient in terms of propellant expenditure, trajectory control is 
enhanced, since the local vertical (local horizontal) may be easily 
sensed. Past the knee of the curve, 0 increases at an average of 
about 3° per day. 

Radial velocity is shown as a function of time in figure 5(b), and 
the similarity to the thrust-angle program is noted. The vehicle's orbit 
is essentially circular up to 100 days at which point the eccentricity 
is only 0.028. A characteristic of the function u(t), which is entirely 
masked out in figure 5(b), is its oscillatory nature over the flat region 
of the curve. A close examination of the digital -computer results, in 
a region where output is called for at every integrating step, shows that 
u oscillates with an increasing period and decreasing amplitude as time 
increases. The period is approximately equal to the orbital period 
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during a given revolution. Initially, the amplitude of oscillation is 
of the order of several meters per second. At t = 25 days, the ampli- 
tude has decayed to a small fraction of 1 meter per second. A corre- 
sponding oscillation is also exhibited by the thrust angle and the 
eccentricity. 

The angular velocity decreases with time by almost four decades, 
as shown by figure 5(c). An approximation of the average orbital period 
may be made^ at any time by 2it/<o. For example, at t = 15 days 1 rev- 
olution is made in about 0.1 day, while at t = 82 days one revolution 
requires about 1 day. Figures 5(d) and (e) show the time histories of 
radial and angular positions, respectively. 

Initial-condition variations . - With reference to equation (6), the 
effect of errors in velocity, position, and mass at t * t Q = 0 on these 
same quantities at tf = 139 days is given by A(t 0 ) . The elements of 
this matrix have been obtained from the adjoint-equation solution and 
are given as 



&u Q 

&co 0 

&r Q 

&9o 

5mo 

5uf 

7.39X10 -4 

7 .20X10® 

1.99X10 -3 

0 

-2.13 

&o>f 

5.69X10” 13 

-2.01K10" 3 

-5. 43X10 “ 13 

0 

5.19X10 -10 

Srf 

1.51X10 2 

1.44X10 15 

4.00X10 3 

0 

-3.80X10® 

5cpf 

-4.41X10 -4 

-1.23X10 7 

-3.41X10 -3 

1 

0.768 


0 

0 

0 

0 

1 


(23) 


Consider, for example, a variation in the initial-circular-orbit 
altitude of 1 10 kilometers, that is, &r Q = 10 4 meters. Assume also that 
5u 0 = &cp 0 = 5mo = 0. Since 



6r 0 = - 2 .08X10“® radian/sec 
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From the second and third columns of A(t 0 ) the following are obtained; 

Buf = 5.00 m/sec 
Boof = -1.25X10”^ radian/sec 
Br-p = 1.00X10 7 m 

Scpj = -8.49 = -2.21 radians = -127° 

The final errors Suf, Sof> and Br^- are each of the order of l/2 percent. 
If there were no error in the angular position, a transformation to 
heliocentric coordinates would show an extremely small perturbation on 
the initial heliocentric conditions. This is not the case, .however, 
since cpp is the parameter that is related to a heliocentric reference 

direction and is in error by -127°. The perturbed escape trajectory 
may be visualized by a 127^* clockwise rotation of the reference trajectory 
in figure 2. Qualitatively, the result is a very significant error in 
heliocentric coordinates. 


The preceding example has served to focus attention on the extreme 
sensitivity of final angular position relative to the sensitivity of the 
other state variables. There is, however, an essential difference in. 
the nature of the sensitivities. That is, the magnitudes of Suf, 5a)f, 
and Br^ increase proportionally with initial-condition errors, whereas 
the principal magnitude of Bcp f varies between 0 and rt in a cyclical, 

triangular fashion. As an example, consider each initial-condition 
error independently and compute the error magnitudes that cause cp f to 

be rt radians out of phase. From the fourth row of the matrix A(t 0 ) 
if Bcpf = ± (2n + l)rt, where n = 0, 1, 2, . . . , 


5u n = + — ^ — j = + (2n + 1)(7.10X10 3 ) m/sec 

4. 41X10 

go) = + — — — =■ = + (2n + l)(2.55X10“ 7 ) radian/ sec 
° 1. 23X10 1 

SCD _r> O 

6r 0 = +■ = T (2n + l) ( 9 . 21X10^) m 

3. 41X10 


&9o = 1 



(2n + l)rt radians 


Bcp-f 

a™ = + , . - L . - ± 
° ° “ 0.768 


(2n + l) (4.09) kg 
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Th6 minimum errors are found when n = 0. From the previous results, 
it may he concluded that &Uq is several orders of magnitude greater 
than any reasonable error that may be expected. In contrast, a 0.04- 
percent error in or a 0.1-percent error in either the initial orbital 

altitude or mo is sufficient, if left uncorrected, to reverse the di- 
rection of the escape asymptote . 

Thrust-vector variations . - Consider next the sensitivity of velocity 
and position components at the nominal escape condition to variations 
in thrust magnitude and direction. This information is expressed by the 
integral terms in equation (6), where the elements of the weighting 
matrix W are given by equations (15) and (16) and are plotted against 
time in figure 6. Each weighting function, except w 41 , exhibits an 

oscillatory characteristic the amplitude and period of which increase 
with time. As in the case of the radial -velocity characteristic, the 
period of oscillation approximately coincides with the orbit rotational 
period of the vehicle. In figure 6, only the last cycle of the weighting 
functions is shown; however, the envelope of oscillation is plotted. 

Qualitatively, the weighting functions illustrate the relative sen- 
sitivity to thrust-vector errors during any two arbitrary t im e intervals. 
Figures 6(b), (d), (f), and (h) show that final velocity and position 
components are relatively insensitive to first-order thrust-angle errors 
during the first half of the escape phase, whereas figures 6(a), (c), and 
(e) show a significant and essentially constant sensitivity to thrust- 
magnitude errors during this same interval. In particular, figure 6(g) 
illustrates a result that could easily have been predicted, namely, that 
9f is most sensitive to SF occurring during the early tightly wound 

spirals and that the sensitivity is essentially zero during the last 
spiral turn (t > 120 days). 

A convenient quantitative measure of the sensitivity to thrust errors 
is obtained by considering SF and Bp to be constant; thus, the time 
integrals of the weighting functions yield the desired sensitivity. The 
results are presented in the following matrix form, where SF and BP 
are measured in newtons and radians, respectively: 



w. 


i 3 


dt 


SF &P 


Sup 

&0D f 

Sr f 

BcPf 


3.73X10 3 

11.7 

-9.11XL0" 7 

-2.77X10" 7 

6 . 68X10® 

7.32X10 7 

-1.35X10 3 

-2.92 


(24) 
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Consider the effects of a thrust -magnitude error first. Since an 
appreciation for the serious consequences of a large error in cp-p has 

been gained, the sequence of 5F values, which causes cp_^, to be it 
radians out of phase, are computed. From equation (24) 

SF a ± ( 2n + - ± (2n + l)(2.33X10" 3 ) newton n = 0,1,2,. . . 

1.35X10 3 

The minimum value of 5F required is thus 2.33X10 -3 newton, or a 0.1- 
percent error. Since it is unlikely that such accurate thrust control 
is possible, it must be concluded that the requirement of guided flight 
is of major importance for low-thrust missions. 

The effects of a ±0. 1-percent error in F on the velocity and 
radial -position components at t-p are 

Buf = ± (3.73X10 3 )(2.33X10“ 5 ) = ±8.70 m/sec 

Bcup = ± ( -9 .11X10" 7 ) ( 2.33X10 -3 ) = +2.12X10" 9 radian/ sec 

Br f = ± ( 6 . 68X10 9 ) ( 2 . 33X10" 3 ) = ±1.56X10 7 m 

Each of these final variations when divided by its respective reference 
quantity represents an error of less than 1 percent^ therefore, the 
resultant perturbation on the heliocentric trajectory is due almost 
entirely to Scp-p = it. 

For discussion of the effects of a thrust-angle error, an error 
of ±10 milliradians, or about l/2°, is assumed as a typical order of 
magnitude. The final velocity and position errors that result are then 
given by equation (24) as 

5Uf = ±0.117 m/sec 
&o^> = +2.77X10 -9 radian/sec 
Br f = ±7 .32X10 5 m 
Scp f = +0.0292 radian = +1.67° 

which are all first-order effects. The error magnitudes due to the first- 
order variation of p appear to be quite small j however, before a con- 
clusion can be drawn the effects of second-order variations of 8 must 
be computed. From equation (22), 



18 


8xJ(tf) = 5xj_(tf) due to 


S0 2 



In this study, F^ eg /2 is of the order of unity. Also, computer results 
show that A i5 /c is small compared with w il for all values of i. 

Thus, the second-order effects are very well approximated by the integrals 
of Wj^L, which are given in the first column of the matrix equation (24), 
multiplied by S(3 2 = 10 -4 : 


Su£ = -0.373 m/sec 
5cq£ = 9.11X10"-*-! radian/sec 
5r^, = -6.68X10 5 m 
Scp^ = 0.135 radian = 7.74° 

A comparison of the first- and second-order effects shows that second- 
order effects are of the same magnitude or greater except in the case 
of angular -velocity errors. The total errors in velocity and radial- 
position components are quite small. The total angular -position error 
is less than 10°. If, however, 5(3 were increased to 2.8°, the final 
angular -position error would increase to about 180°. 


Trajectory Sensitivity - Heliocentric -Transfer Phase 

R eference trajectory . - The time required to transfer between plan- 
etary orbits was chosen as 223 days, and the trajectory was optimized 
to minimize the propellant expenditure. The constraint imposed here 
was that the thrust must be equal either to the design value or to zero. 
Characteristic variables of the reference trajectory are plotted in fig- 
ure 7 where the time scale chosen has been reinitialized to zero. The 
optimum- thrust -angle program is shown in figure 7(a). Thrust direction 
is outward with an average rate of change of 0.86° per day during the 
first powered period and inward with an average rate of change of 0.56° 
per day during the second powered period. Velocity and position time 
histories are shown in figures 7(b) to (e). 


Initial-condition variations . - Digital-computer solutions of the 
adjoint equations (8) to (12), which correspond to the boundary conditions 
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A(tf) = unit matrix, are plotted as functions of time in figure 8. The 
effects of initial-condition variations on final conditions are given 
"by - A(t 0 ) : 



Su 0 

Bcd 0 

Sr 0 

Scpo 

Sm 0 

Su.^» 

-0.366 

2.63X10 11 

6.46X10 -7 

0 

-0.149 

&0>£ 

-1. 84X10 ~ 12 

-3.48 

-8.07X10- 18 

0 

8.07X10- 12 

A(t 0 ) = Brf 

2.2CX10 6 

4.57X10 18 

10.9 

0 

-1.2 7X10 7 

S9 f 

-1. 02X10 “ 4 

-3.11X10 7 

-8.44X10- 11 

1 

. 

1.00X10" 4 

Bm^ 

0 

0 

0 

0 

1 


(25) 


Variations in initial velocity, position, and mass will be due to errors 
incurred during the escape-spiral phase. Expressions are derived in 
appendix B that relate the two sets of errors. Suppose, for example, 
at the termination of the escape phase that the angular position is in 
error by only 10 milliradians, while all other errors are zero. With 
the use of equations (B6), (B7), (B12), and (B14) , the initial- 
heliocentric -velocity and -position errors are found to be 

5u 0 = -15 . 8 m/ sec 

Su>o = 2 .47X10 radian/sec 

Sr Q = -1.78X10 7 m 

Scp 0 = 4. 50X10 radian 

&m Q = 0 

Then from equation (25) the final -condition errors will be 

6U|> = 0.770 m/sec 

5a>f = 8. 69X10”^ radian/sec 
5r f = -1.16X10 8 m 
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Sq)|> = 2.39X10“^ radian; i^cp-p = 5.43X10® m 


5mf = 0 

With reference to figure 3, the resultant errors show that at the nominal 
final time the vehicle is both below and ahead of its nominal position, 
and its motion relative to Mars is essentially parallel to the nominal 
capture trajectory. Hence, the incoming asymptotic displacement is in 

error by some 1.16X10® meters. The requirement for midcourse corrective 
guidance to compensate for the escape errors is evident. 


Thrust-vector variations . - The elements of the matrix W are 
plotted as functions of time in figure 9. Recall that the nominal coast- 
ing period extends over the interval 27.7 < t < 187.7 days. Since it 
is reasonable to assume that no thrust errors will occur during this 
interval, the weighting functions are of interest here only during the 
nominal powered periods. The significant result illustrated by figure 9 
is that the final velocity and position components are relatively more 
sensitive to thrust-vector errors that occur during the first powered 
period. In particular, the radial -position perturbation due to either 
5F or SP is relatively 15 times greater during the first powered 
period, while the angular -position perturbation due to SF is relatively 
25 times greater. 


The time integrals of the weighting functions over the powered 
regions of flight have been evaluated as 


/ 

•'Tin 


w ij dt 


powered 

periods 



SF 

S3 

5u.f 

211 

-604 


-1. 14X10-® 

2. 45X10-® 

Sr^ 

1.80X10 10 

-2.15X10 10 


-0 .142 

0.0384 


(26) 


As an example, the velocity and position perturbations resulting from a 
1-percent error in thrust magnitude or a 10-milliradian error in thrust 
angle are 



&F = 2.52X10"® newton 5(3 = 10~® radian 


6 Up = 4.90 m/sec 


6up = -6.04 m/sec 

5to f = -2.64X10" 10 

radian/ sec 

= 2.45X10 - -^ radian/sec 

5r f = 4.17X10® m 


5r f = -2.15X10® m 

5cp f = -3.29X10" 3 

radian 

bqxp = 0.384X10 -3 radian 


The resultant errors and, in particular, the position errors are 
of sufficient magnitude to require midcourse guidance corrections . 


Trajectory Sensitivity - Capture -Spiral Phase 

Reference trajectory . - Variables that describe the capture-spiral 
trajectory are plotted as a function of time in figure 10. The thrust 
angle (B, corresponding to the energy-decreasing tangential-thrust 
program, is shown in figure 10(a) in which p is shown as a second 
quadrant angle, since the thrust vector in this case is parallel but 
opposite to velocity vector. As in the case of the escape-spiral 
trajectory, the thrust direction is essentially circumferential over a 
large region. In this region, t > 25 days, the vehicle travels in a 
nearly circular orbit with a decreasing semimajor axis. These results 
may readily be inferred from the radial -velocity and -position curves 
shown in figures 10(b) and (d) . Again, as in the case of the escape 
spiral, a characteristic of the function u(t), which is entirely masked 
out in figure 10(b), is its oscillatory nature over the flat region of 
the curve. The period and amplitude of oscillation decrease and increase, 
respectively, with time, although the largest amplitude is quite small. 

Variational parameters of final satellite orbit . - It may be more 
illuminating to consider the sensitivity of the final satellite orbit 
about Mars in terms of orbital parameters such as semilatus rectum p, 
semimajor axis a, eccentricity e, perigee rp, or apogee r a - If the 

orbital orientation is not of interest-, any two of these parameters will 
define the orbit. The parameters a f and p f may be expressed in 

terms of the radial -position and -velocity components as 

a f = — 

2 u f + (lyx ^) 2 
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r fCo| 

Pf = s%- 

The eccentricity, perigee, and apogee are given as 

e^ a l - — 

f a- 


— ap(l - ep) 


r a,f = a ^ 1 + e f) 

When the reference orhit is circular, the following special relations 
occur: 

Uf = e f = 0 


a.f - Pf - - r a ^£ r f 


V f r f “ f - l rf 


Expansion of ap and pp in a Taylor series keeping only first- and 

second-order terms and utilization of the previous special relations 
yields 


lap — 


f k „2 


+ + /20 


+ «*p + 77 K + 72" M + 77 5r f + f 


Ap^ = 


fk 2 , / 6 V 2 / 8 


Variations in eccentricity, perigee, and apogee may then he calculated 
from 

(29) 


Arp p = Aap - (rp Aep + Aap Aep) 
Ar a ^p = Aap + (rp Aep + Aap Aep) 
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Note that, even for very small variations in velocity and position, equa- 
tions (29) to (31) are essentially nonlinear when expressed in terms of 
these variations. This result is due to the fact that the final reference 
orbit is circular. 

Initial-condition variations . - Solution of the adjoint equations 
for the capture-spiral trajectory yields the following A(t 0 ) matrix: 


SUf 


SoDf 


A(t 0 ) = 8r f 


8<P f 


5mp 

(32) 

The magnitudes of the sensitivity coefficients in equation (32) indicate 
the extreme sensitivity of the capture trajectory to errors in initial 
conditions. A sample calculation illustrates -this point. Consider a 
situation where the vehicle is at the influence sphere (&r c = O'), hut 

the velocity components are in error hy 5u 0 = 1 meter per second and 
5<n 0 = -10 -9 radian per second. This corresponds to the initial velocity 

vector having an error of about 1 meter per second in magnitude and 1 
milliradian in direction. Assume also that Sm Q = 0. From equation (32), 

the final variables will be in error by 

5u f = -24.2 - 997 = -991 m/sec 
Boop = -2 .02X10” 4 - 2 .01X10 = -4.03X10 -4 radian/sec 

5r f = 5.73X10 5 + 5.77X10 5 = 1.15X10 6 m = 1150 Inn. 

These errors represent a significant perturbation of the trajectory. 
Evaluation of equations (27) to (30) shows the eccentricity to be about 
0.5 and the perigee to be well within the Martian surface. Actually, 
the errors are large enough to invalidate the linear perturbation approach. 

An interesting result may be found from the matrix A(t 0 ) if in 
any given column the element of the second row is divided by the element 


5u 0 5<n 0 5r 0 &<Po ° m o 


-24.2 

9.67X10 11 

5.37X10" 5 

0 

7.16 

-2.02X10' 4 

2 .01X10 5 

-l.ioxio- 10 

0 

4. 87X10 “ 5 

5.73X10 5 

-5.77X10 14 

0.311 

0 

-1.39X10 5 




1 






0 

0 

0 

0 

1 
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of the third row. These ratios are very nearly the same, about 

-3.5X10 -10 . Within the validity of the first-order approximation, this 
result means that 5cU|. is not independent of Sr^; that is, no matter 

what set„©f initial errors leads to a given 5r f , Sco f « -3.5X10 -10 &r f . 
This common ratio is not unique to the time t Q = 0, since an examination 
of the matrix A(t) shows approximately the same result over a large 
portion of the trajectory. 

A further result is that the quantity 2co^/r f is equal to 
3 .54XQ.0-3-0; thus, Sco^ -{Za^/xf )&rp . This result is not explained in a 
rigorous manner; however, a relation of this form and order of magnitude 
may be deduced from the following argument. If the errors 5u f , &CD f , 
and Sr f are quite small, the final perturbed orbit is very nearly 
circular. Since for a circular orbit 


CD-p 



the relation between Soi|= and Sr^. required to maintain a circular 
orbit of slightly different size is 

*■* - -I 


The 25-percent discrepancy in the coefficient &r f may be ascribed to 
the fact that the final perturbed orbit is not exactly circular. 

Thrust-vector variations . - The sensitivities of final -velocity 
and -radial-position components to variations in thrust magnitude and 
direction are given by the weighting functions, which are plotted in 
figure 11. As in the case of the escape-spiral .phase, the functions 
exhibit an oscillatory characteristic, the envelope of which is shown. 

The major qualitative result to be gained from figure 11 is one that 
may have been anticipated: the maximum sensitivity to thrust errors 

occurs in the vicinity of t G when the vehicle is far from the planet 

and moving asymptotically, and the sensitivity decreases to a negligible 
amount as the final satellite orbit is approached. The guidance impli- 
cations are clearly understood. Trajectory perturbations due to initial - 
condition errors are most efficiently corrected early in the capture 
phase, provided, of course, that accurate guidance information is 
available to the vehicle at such large distances from the target planet. 
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A quantitative measure of the sensitivity to thrust variations is 
provided hy the time integrals of the weighting functions : 



Sr-p 


5F 5p 


-9.03X10 3 

-3.44X10 5 

-6 .15X10- 2 

-7.83X10-2 

1.75X10® 

2 . 23X10® 


(33) 


A simple calcination shows that a significant perturbation of the final 
satellite orbit results from thrust errors as small as 0.1 percent in 
magnitude and 1 milliradian in direction. For SF = 2.32X^.0~^, the 
final -velocity and -position errors are -20.9 meters per second, 

-1.43X10" 4 radian per second, and 4.05X10 5 meters, respectively. For 
5(3 = 10 - ®, the final errors are -344 meters per second, -0.783X10 -4 ra- 
dian per second, and 2.23X10 3 meters, respectively. In the last case 
the eccentricity of the perturbed orbit is about 0.1, and the perigee 
altitude has been decreased from the nominal 927 to 390 kilometers. 


For reasons previously discussed, the second-order effects of 

2 

thrust-angle variations may be approximated by the product of Sp and 
the integrals of wqq given in the first column of equation (33). 

Since the integrals of wpi are of smaller magnitude than the integrals 
of w i2 , the second order effects are negligible for Sp <0.1 radian. 

Note, from equation (33) that the ratios of the elements taken from 

the second and third rows are both equal to -3.51X10 Thus, f rom the 

arguments of the previous section So-p ~ -(2a)|*/r-f )5r.p for any small 

values of constant SF and 8P • 


Examples of Simple Guidance 

Escape guidance in one variable . - The basic result of the sensitiv 
ity analysis of the escape -spiral, trajectory is that a very small pertur 
bation in thrust magnitude is sufficient to rotate the asymptotic escape 
direction through a large angle. Specifically, if a 0.1-percent system- 
atic error exists and no corrective maneuver is made, the escape 
asymptote will lie on the sun’s side of the Earth's orbit and will be 
directed, opposite to Earth's orbital velocity. This condition is 
equivalent to a 25-percent loss in heliocentric energy, which would most 
likely result in a mission failure. The objective of escape guidance 
is to prevent such a condition from occurring. 
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As an example, consider a simple guidance scheme that employs 
constant thrust AF as the control variable and a single -variable guid- 
ance criterion, namely, to null A constant-thrust-magnitude error 

(unknown to the vehicle) of 0.1 percent is arbitrarily assumed as the 
uncontrolled perturbing function, and repetitive corrective action is to 
be taken at 20-day intervals based on perfect trajectory determination. 
The calculation is simplified by restricting the control interval At 
to within 1 day, where the thrust weighting functions are essentially 
constant over this interval. The details of the thrust -control system 
are left unspecified, and results are given in the form of thrust 
impulse AF At measured in newton-seconds. Depending on the control 
required, then, consideration may be given to the possibilities of throt- 
tling the low-thrust engine, shutting it off completely, or using a 
medium-thrust chemical rocket. 

The guidance equation involving angular position only may be ex- 
tracted from equation ( 6 ) . At the beginning of any control interval 
t n the final-angular-position error predicted from measurements is 
denoted as 8cp-f(t n ) . The result of equation (2l) is substituted for 

the integrand term in equation (6)j hence, the guidance equation may be 
written as 


Arpf (t n t At n ) 


Ap>f ( t n ) 


+ 



W4qAF n dt 


= Acp f (t n ) + w 41 (t n )AF n At n 


Since Ap f (t n + At n ) = 0 


is desired, the control thrust impulse is 


^n ^n " 


Ap f (t n ) 


If corrective action were taken at the previous control instant t n _q 
null Apf, Atp|>(t n ) is due to the perturbing function 5F acting over 
the interval (t n _^, t n ) . For the purpose of the calculation herein, 


AF n Atn — 





to 


Figure 12(a) illustrates the results of this simple -final -value 
guidance scheme. A constant-thrust perturbation 5F = 2. 32X10-3 newton 
(0.1 percent) and five control intervals at 20, 40, 60, 80, and 100 days 
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are assumed. The uncorrected Acp-p characteristic is shown for compar- 
ison. Although the last corrective action was taken at 100 days, the 
angular-position error increases to an acceptable final value of less 
than 1° during the interval 100 to 139 days. Since a positive pertur- 
bation was assumed, the control impulse must be negative to make up for 
the lag in angular position. This impulse increases in magnitude 
slightly with each successive correction from 5.08X10 3 to 8.34X10 3 newton- 
seconds. Consider the correction at t = 20 days. If the low-thrust 
engine is shut off completely (AF = -2.32), the shutoff time At is 
2 .18X10^ seconds, or about 0.6 hour. If the engine is throttled down 
by 5 percent, the control interval is about 12 hours. 


The effect of the guidance maneuvers on the uncontrolled final 
variables SUf, Sco^, and Sr^ is shown in figures 12(b) to (d). Control 

stability results in each case, and a comparison with the uncorrected 
characteristics shows a significant reduction in the final errors . 

Midcourse guidance in two variables . - As a second example of how 
the fundamental guidance equation may be used to prescribe corrective 
maneuvers, guidance action taken during the nominal coast period of the 
heliocentric-transfer phase is considered. The coast period begins 
27.7 days after heliocentric injection, and this time instant is denoted 
t c . If the powered period is extended over the interval (t c , t c + At) , 
where At is arbitrary but less than several days, and if the thrust - 
angle correction A(3 is constant during this 'interval, any two of the 
final -velocity and -position error components may be nulled. The 
position errors were chosen to be nulled. 


The guidance equations involving position errors only are taken 
from equation (6) and are written as 


5r f (t c + At) = 6r f (t e ) + 


/ t 0 +At 


f 3 dt 


scp f (t c + At) a s<p f (t c ) + / f 4 dt 


/ t 0 +At 


where 5r-f(t c ) and &<Pf(t c ) are the final -position errors predicted at 

time t from navigational measurements, and the integrands f„ and 
c o 

f 4 are given by equation (20). For the purpose of calculation, the 
integrands may be approximated by dropping the terms ^3.5/c, since 
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results have shown them to he quite small in comparison with the other 
terms. Therefore,, since &r f (t c + At) = 5cp f (t c + At) = 0 is desired, 

and At is small. 


&r f (t c ) + P, 


des 


At 


w 32 (t c ) 

|W 3 i(t c ) co s AP + ^ sin AP 


F 


= 0 


and 


Scp f (t c ) + F 


des 


At 


w 


(t )cos AP + 


w 42( t c) 


41 ' c 


F 


sin AP 


= 0 


The solutions for AP and At are 

w 4l(' t c )S r f(' fc c J " 


tan AP = 


w 


F S®f(h c ) F 


42 ( t c ) 


5r f (t c ) 


At = 


v- 


-5r f (t c ) 


i , * w 32 (t ) 

^des | w 31^c^ cos ^ + f sin AP 


J 


where the quadrant of Ap is chosen such that At is positive. The 
effect of the corrective maneuver on the final -velocity errors may he 
determined from 


6 uf(t c + At) = &Uf(t c ) + Fq es At 
5co f (t c + At) = Soj f (t c ) + %es At 


|Wqq(t c ) co S AP + 

|^21 (he) 003 + 


*12^ 

F 

w 22 (t c ) 

F 


sin Ap 
sin Ap 


As a numerical example, assume the final errors determined at time 
t c to he due to the residual escape-guidance errors and to a constant 
0.1-percent-thrust-magnitude perturbation SF acting over the interval 
t 0 = 0 to t Q = 27.7 days. The final-escape errors from the previous 

section are transformed to initial heliocentric errors with the use of 
equations (B6), (B7), (B12), and (B14) . Then, from the results of the 
heliocentric -trajectory-sensitivity analysis 


5uf(t c ) = 0.240 m/sec 


5oif(t c ) = -1 .58X10"-'-® radian/sec 
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Sr f (t c ) = 2.08X10 8 m 
5cpp(t c ) = -3.93X10"^ radian 
The midcourse corrective maneuver requires that 

A(B = 3.91 radians = 224° 

At = 5.66X10 4 sec = 0.656 day 

The actual thrust direction is found hy adding Af3 to the reference 
thrust direction at t G (see fig. 7(a)); hence , 

P = 0.850 + 3.91 = 4.76 radians = 273° 

The mass loss due to the corrective maneuver is 

Am = m At 

= (-5.68 kg/day) (0.656 day) 

= -3.73 kg 

The final -velocity errors after the correction "become 

5Uf(t c + At) = 9.83 m/sec 
&ocif(t G + At) = 2.25X10"-^- radian/sec 


CONCLUDING REMARKS 

A typical low-thrust Mars orbiter mission using the Snap-8 power- 
generating system is studied from a guidance viewpoint . The mission 
trajectory is divided into the three characteristic phases, namely, the 
outward escape spiral, the heliocentric transfer, and the inward capture 
spiral. For each trajectory phase, the sensitivity of the final trajec- 
tory state (velocity, position, and mass) due to errors in the initial 
trajectory state and to thrust-vector errors is determined. This infor- 
mation is expressed by the fundamental guidance equation, the derivation 
and solution of which is based on linear perturbation theory and the 
method of adjoint functions. In addition to providing a means of inves- 
tigating the perturbative effect due to a number of error sources, which 
is the major objective of this report, the guidance equation may also 
be used to determine requirements for corrective guidance. 

The escape-spiral trajectory is designed so that the hyperbolic 
escape asymptote is pointed in a prescribed direction relative to the 
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Earth-sun line. In other words, the last spiral turn must he oriented 
properly. Since the vehicle makes several hundred revolutions about the 
Earth in the process of escaping, it might be expected that the escape 
direction would be highly sensitive to errors. This expectation is amply 
verified by the results of the sensitivity analysis. For example, a 0.1- 
percent error in either the initial-orbital altitude or the initial mass, 
if left uncorrected, is sufficient to reverse the direction of escape. 

A minimum error of 0.1 percent in the thrust magnitude acting over the 
entire trajectory will also result in a 180° misorientation. The effects 
of thrust-angle errors during the escape phase are found to be nonlinear 
for errors larger than a fraction of 1°, so that second-order terms are 
required. Results show, for example, that thrust-angle errors of l/2° 
and 2.8° acting over the entire trajectory cause the escape direction 
to be in error by about 10° and 180°, respectively. Because of the 
small control errors involved, specifically the thrust -magnitude error, 
the vehicle must have the capability of corrective-guidance programming. 

The sensitivity analysis of the heliocentric transfer phase shows 
that the nature and magnitude of trajectory perturbations are not unlike 
those for free-fall trajectories. The exception, of course, is that an 
error in the vehicle's mass will perturb the low-thrust trajectory. An 
error of 10 kilograms at the initial -trajectory state results in a final- 
position error of several hundred thousand kilometers. Either a l/2° 
error in the escape direction, a 1-percent error in thrust magnitude, or 
a 1/2° error in thrust angle, will result in final position and velocity 
errors of the order of several hundred thousand kilometers and 10 meters 
per second, respectively. These errors are of sufficient magnitude to 
require midcourse corrective maneuvers . 

The capture spiral, like the escape spiral, is highly sensitive to 
errors both in the initial trajectory state and in thrust-vector control. 
These two phases of the mission are essentially duals of each other. In 
the escape phase, the final velocity and radial position are not too sen- 
sitive to guidance errors. The important parameter is the final angular 
position, which is strongly related to the escape direction. In the cap- 
ture phase, however, the final angular position is of little consequence 
compared with the size of the final satellite orbit . Hence , the final 
velocity and radial position are the important parameters, and results 
have shown these to be very sensitive. For example, if the initial- 
velocity magnitude and direction were in error by only 1 meter per second 
and 1 milliradian, respectively, and if the reference thrust program 
were followed exactly, the vehicle would spiral down into the Martian 
surface instead of establishing the nominal 927 -kilometer circular orbit. 
A constant -thrust -angle error of only 1 milliradian acts to decrease 
the final perigee altitude by about 400 kilometers . 

In summary, the results of this analysis indicate that the low- 
thrust trajectory is extremely sensitive to relatively small error 
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magnitudes. It should not be concluded, however , that accurate guidance 
is unachievable. As the trajectory is affected by error perturbations, 
it is likewise affected by controlled perturbations, that is, by correc- 
tive thrust programming. The proper conclusion to be drawn from this 
study is that the electrically propelled space vehicle incurs the job 
of repetitive trajectory determination and corrective guidance maneuvers. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 16, 1962 
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APPENDIX A 


PERTURBATION ANALYSIS 

The motion of the vehicle during each phase of the mission is de- 
scribed by a set of nonlinear differential equations. Consider that a 
solution of such a set corresponding to a specified thrust program and 
satisfying prescribed boundary conditions has been obtained and is 
termed the reference solution. If the major assumption is made that 
extraneous perturbations (e.g., guidance errors) are sufficiently small 
so that the actual vehicle trajectory does not vary significantly from 
the reference trajectory, it is possible to study these variations and 
the required corrective maneuvers by linear perturbation techniques. 

The perturbed differential equations of motion represent a linear system 
with time-varying coefficients, which are expressed as known functions 
of position, velocity, thrust and gravity forces, and mass along the 
reference trajectory. 

The following analysis shows the derivation of the linearized 
system equations and the use of adjoint methods in obtaining the varia- 
tional solution and expressing the fundamental guidance equation. In 
the following discussion, the symbol 5 is used to represent small 
variations from reference quantities, and matrix notation is used for 
the purpose of compactness and ease in algebraic manipulation. 


System Equations 


During each phase of the mission, the vehicle is assumed to travel 
in a vacuum under the influence of an inverse -square central gravita- 
tional field in addition to its own thrust acceleration. A two- 
dimensional geometry is used, and the vehicle motion is described in a 
rotating polar coordinate system centered at the appropriate centr al 
body. With reference to figure 1, the differential equations of motion 


the flight path are given as 


2 GM 
u = r<xr — 

r^ 

F 

+ — sin (B 
m p 

(1) 

. -2uoo . 

CD = + 

r 

F 

— cos B 
rm 

(2) 

r = 

u 

(3) 

cp = 

00 

(4) 

m = - 

. I 

c 

(5) 
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Differentiation with respect to the independent variable time is denoted 
by a superscribed dot. The state variables of interest are radial 
velocity u, angular velocity in, radial position r, angular position 
9 , and mass m. The propulsion or control variables are the thrust 
magnitude F, and the thrust angle (3, which is measured with respect 
to the local horizontal (circumferential) direction. The effective jet 
velocity c is equal to the product of specific impulse and a conversion 
factor, and is considered constant herein. 

At this point, it is convenient to define the following column 
matrices or vectors: 



and to write equations (l) to (5) as 

x ± = g i (x 1 , • • X 5 , y 1 , y 2 ) i = 1, • • ., 5 (A2) 

The perturbed system of equations is obtained by taking the first-order 
variations of (A2) 



j=l 


which is, in matrix notation, 

Sx - A&x = B&y (A4) 

where A and B are (5 by 5) and (5 by 2) coefficient matrices, respec- 
tively, whose ij th elements are given by dg^/dx^ and 

respectively) thus. 
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Note that A and B are time-varying matrices evaluated along the 
known reference trajectory. 


Solution by Adjoint Methods 

The solution of equation (A4) can best be facilitated by the method 
of adjoint functions, as suggested in reference 3, by introduction of a 
5 by 1 vector of Lagrangian multipliers Aj_, which satisfies a system 

of equations defined to be adjoint to equation (A4) : 

\ + A^Ai = 0 (A7) 

where 


A± = 


In equation (A7,) A is the transpose of the matrix A. The desired 
relation between equations (A4) and (A7) is obtained by taking the __ 
scalar product of (A4) with A-j_, the scalar product of (A7) with 5x, 
and adding the results: 
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Af * (^jT 8X^ - \ • (ABx) + A-^ * 5 x + (^A*) • Sx = A ± • (B5y) 

It is easily shown, however, that 

\ • (A5x) = ^A-j^ • 5x 

Aq • (B6y) = (b^A^ * &y 

ik (X i • 63S) = A ' (It w ) + (it ^ ' 635 

Therefore 

JjT (Ai * Sx) = * 5y (A8) 


Define a 2 hy 1 weighting vector wq : 


W-t 


Wii 

_w i2 _ 


= 


(A9) 


and integrate both sides of equation (A8) between the general time t 
and the nominal final time t^ i 


Ai(tf) * 5x(t|*) = Aq(t) • 5x(t) + 




• &y)dt 


(A10) 


Equation (A10) expresses a linear function of variations in final 
conditions in terms of variations in state variables at any point along 
the path and the integrated effect of future thrust variations over the 
path. For the purpose of analysis it is desirable to separate the effect 
of variations on each quantity individually. This is possible by a 
proper interpretation of the boundary conditions on the Lagrangian multi- 
pliers. For example, let the variation Sx^ be of interest at the final 

time t f . If the following boundary conditions at t f are specified, 

Aii(tf) = 1; Aij(tf) =0 for j ^ i (All) 

the adjoint equations (A7) can be integrated backward from tf to yield 
numerical values for the multipliers at each time instant t. The varia- 
tion of interest is then given explicitly from (A10) : 
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Sxi(tf) = Ai(t) • 5x(t) + / (wi • By)dt 

't 


r 


(A12) 


If this procedure is followed for i = 1, . . .,5, the complete solution 
for the variations in final-state variables may he expressed in the 
matrix form 


5x(t ;£ .) = A(t)Sx(t) + 


r 


(W5y)dt 


(A13) 


where the rows of A and W are made up of the elements of A and 
w ± , respectively: 



V 


‘ A ll 


• 


A 21 

A = 

• 

— 

• 


_v 


_ A 51 


A 


‘15 


A, 


25 


A 


55 


; W 


¥ n 


Wc 


W' 


'll w 12 


[_ w 51 W 52J 


(7) 


Because of the form of equation (7), the elements of A can he inter- 
preted as partial derivatives or sensitivity coefficients; thus. 


Aij(t) = 


dx ± (t f ) 

dXj(t) 


ji — Ij • • •> 
y = 1 j • • • ) 


(All) 


Wow that the general results have been established with the help of 
matrix notation, it would be well to write out equations (A7) and (A9) 
in terms of the variables of interest in this study. The set of adjoint 
equations is determined from (A7) and (A5): 


A il = (|r) A i2 " 7*13 
A i2 = - (2rco)A il+ (z fj 


A i2 " A i4 


^ / 2GM 2K 

\3 = - + “ ) A il 


' 2un> F cos 3 
ir^ r^m 


A i2 


A i4 = 0 


A i5 - 


F sin [3 

? 
nr 


*11 + 


( 8 ) 

(9) 

(10) 

( 11 ) 

( 12 ) 
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The weighting functions are found from (A9) and (A6): 

w il = 


A il 


m 


. . . A i2 A i5 

sm B + cos B - 

rm c 


■L’ l-\ r, i2 . 0 

w. 0 = — (A., cos p sin p 

i2 mill p r 


(15) 

(16) 


Nonlinear Thrust-Vector Perturbations 


The previous results were derived with consideration given only to 
first-order variations in the state variables and the thrust vector. 

It is possible that relatively large variations in thrust magnitude and 
direction will occur in such a way that the deviation between actual 
and reference trajectories still remains small; that is, the lineariza- 
tion of state variables is still valid. The integrand in equation (A13) 
must be modified from its linear form to account for this situation. 


When exact variations in thrust magnitude and direction (A F and 
Ap) are considered, the right side of equation (A4) becomes 


BSy 


F + AF / n F 

— sin (p + Ap) - — sin p 

m m 

- F + ^ cos (p + Ap) - — cos p 
rm rm 

0 

0 


(F + AF) F 
C + C 


Correspondingly, from the definition of equation (A9), the integrand 
term of equation (A12) becomes 


A . , 

(wj. • &y) -* f ± = ” [ (F + AF) sin (p + Ap) - F sin 0] 
A- p A.? k 

+ — [ (F + AF) cos (p + AP) - F cos p] - — 
rm c 


[ (F + AF) - F] 
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Equation (A13) then becomes 


5x(t f ) = A(t)6x(t) + / (fj:)dt 

•'t 
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APPENDIX B 

ERROR TRANSFORMATION BETWEEN ESCAPE AND 
HELIOCENTRIC MISSION PHASES 

The relations between final escape errors and initial heliocentric 
errors are presented. The Earth-referenced velocity and position com- 
ponents at the termination of the escape phase are denoted as u^p, 
ojfij and cp.Q_ and the sun-referenced components at the initiation 

of the heliocentric phase as uq 2 > t 02 > an< ^ c f ) 02 • Only first- 

order variations in > ^fi^ an< ^ r fl are considered; Jaowever^ large 
variations in cpfi are allowed. 


Position Errors 

The geometric configuration of the Earth-centered reference position 
r fl and the actual position r fl + 6r fl is shown in the following 

sketch. The coordinate frame has i r as a unit vector away from the 
sun,, and i^ as a unit vector in the direction of Earth's orbital 

velocity. 
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The quantities bp and Sp are components of the reference-position 

vector. Since it has been arbitrarily assumed that the hyperbolic 
escape direction is parallel to i^, the asymptotic displacement of the 

escape hyperbola is represented by b-j_ . 

The heliocentric position of the vehicle at time tfp is 

?02 * + 7 fl 

so that 


Sr 0 2 = Sr fl = (&r Q 2 )i r + ( r 02 Scp 02 ).I^ 


The components of 

SrQ 2 , which may easily be found from the sketch, are 

given by 




^ r 02 = 

[ 

r fl L 

r fl (cos - l) + cos Ap fl 5r fl J - 

si sin Apn (Bl) 

r 02 &c P02 

S 1 

" r fl 

[rfl (cos - l) + cos ^Pfp& r fpj 

+ bi sin Apn 




(B2) 

where 


/ r fl“fl\ 

111 = tn [^ ) 

(B3) 



h - hifc) 

(B4) 



r 02 ~ r E + ^l 

(B5) 
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Positive Apfi is measured counterclockwise, so Ap-^ is a negative 
angle, as indicated by the sketch. 


Velocity Errors 

The heliocentric velocity of the vehicle at time t fl is 

^02 = + ^fl 


so that 


where 


SV 02 ' SV fl - feVb + 5(r 02 ffi 02 )l 9 


8 ^ r 0Z°W “ r 02 6 “02 + “02 ST 02 


The magnitude of &Vf]_ may 8e found as follows: 

v fl = u fl + ( r fl“fl^ 

/ufA / r fl ffi fl\ _ , /^fl“fl^ 

&V^ = I — ) Siv. + (— — ) s “fi + 




v fl ~ [VfiJ " fl 


'fl 


V. 


5r 


fl 


fl 


or from (B3) and (B4) 


> 


/&! 


OT fi - fed Su a + + Ivjyy 


&r 


fl 


(B8) 


The angle between the velocity vector and the Earth-centered local 
horizontal is denoted by y in the sketch. The variation in this angle 
may be expressed in terms of 8uj]_, &<&□_, and Si’fi as follows: 


tan y = 


u fl 

r fl“fl 


1 ' 5u fl - ( ■ Uf - ^- )5 (Of! 


^fl^fl, 


1 r fl m fl 


a fl 

^fl^fl 


5r fl 


5y = cos^y 
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since 


cos X = 


r n“fi 

r ti 


&r = 




V' 


Su. 


fl 


fl 


( u fl r fl) 


fl 


&(JD. 


'fl 


i v li ) 


&r fl 

(B9) 

&v f l 

is in part 


due to &T and in part due to Specifically, the angle "between 

the latter vector and i^ is 5T - ^Pfi measured counterclockwise 
positive from + 5V_^ . how the heliocentric components of 5V^ 

may he expressed as 


&Uq 2 — (Vfi + 5Vfi) sin(6Y ~ 

&<r 02 V - (v fi + w fl )=os(8r - - V 


fl 


Since 6 Vfi « V fl and 5 t is very small, the following approximations 
are valid: 

&UQ 2 = V pphT cos /xp fi - V fi sin £cpfi 

&( r o2^2^ = Vfi (cos Apfi - l) + &Vfi cos Zlpff + VfiSy sin Apfq 

With the use of equations (B3), (B4), (B8), and (B9), the previous 
expressions become 


Su 02 “ 


U 1 

, r fl> 


&Uf i ~ ( Sf ) 5c0fi — 


fl^l' 

V r fl > 


Sr. 


fl 


cos APfi - Vfi sin Apfi 

(BIO) 


& ( r 02 CD 02) = < S 1 cos *Pfl + ^1 sin ^Pfl) &u fl 

+ (bi cos Apfi - Si sin 
“fl 

+ (bi cos /^pfi - si sin Z^fi)&rfi + Vfi(cos Apfi 

r fl ' " ~ 


- 1) 


(Bll) 


43 


IT Ap.Q. is very small Ap.Q_ -* 5<p ^ , equations (BIO) and (Bll) 

become 


( bn \ /siO-piX 

SUfi - (s 1 )5(D fl - ( J Sr f] _ - (V fl )8cp fl (B12) 

&( r 02 a ^02^ = ^fl (B13) 

The variation Soc^g is found from (B6) and (B8) 

5a ^2 = L l ■) su fi " f^-Wi 

rfl 02/ \ 02/ 

h-i /avvosA 

+ r fi r 02 ^ 02 " “fl^ 5r fl + ^ Tq2 J BCp fl ^ B14 ^ 
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Figure 1. - Overall interplanetary transfer between Earth and 





Figure 2. - Escape-spiral trajectory. Hyperbolic-escape conditions: time, 139 

days; radial position, 1.899X10 9 meters (298 Earth radii); relative velocity, 
1577 meters per second. 



Asymptote of 
capture hyperbola 


Figure 3. - Capture-spiral trajectory. Hyperbolic-capture conditions: 
radial, position, 1.382x109 meters (418 Mars radii); relative velocity 
1454 meters per second. 
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Figure 4 . - Coordinate reference schematic 
diagram. 
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(a) Thrust-angle program. 
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Time, t, days 

(Is) Radial -velocity time history. 

Figure 5. - Continued. Characteristics of escape-spiral reference 
trajectory. 



Time, t, days 

(c) Angular-velocity time history. 

Figure 5. - Continued. Characteristics of escape-spiral reference 
trajectory. 
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(d) Radial-position time history. 


Figure 5. - Continued. Characteristics of escape-spiral reference trajectory. 


Time, t, days 

(e) Angular-position time history. 

Figure 5. - Concluded. Characteristics of escape-spiral reference 
trajectory. 
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Envelope 


(e) Sensitivity of final radial position to thrust 
magnitude variations. 


Envelope 


Time, t, days 

(f) Sensitivity of final radial position to thrust- 
angle variations. 

Figure 6. - Continued. Thrust-vector weighting functions 
for escape-spiral phase. 
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(h) Sensitivity of final angular position to thrust- 
angle variations. 


Figure 6. - Concluded. Thrust-vector weighting functions for 
escape-spiral phase. 
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(t>) Radial-velocity time history. 

Figure 7. - Characteristics of heliocentric reference 
trajectory. 
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(e) Angular-position time history. 
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Concluded. Characteristics of heliocentric 
trajectory. 
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(fa) Sensitivity of final angular velocity to variations 
in velocity, position, and mass. 

Figure 8. - Continued. Adjoint solution for 
heliocentric-transfer phase. 
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(c) Sensitivity of final radial position to variations 
in velocity, position, and mass. 

Figure 8. - Continued. Adjoint solution for 
heliocentric -transfer phase. 
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(d) Sensitivity of final angular position to variations 
in velocity, position, and mass. 

Figure 8. - Concluded. Adjoint solution for heliocentric 
transfer phase. 







Time , t , days 

(a) Sensitivity of final radial velocity to thrust- 
vector variations. 

Figure 9. - Thrust-vector weighting functions for 
heliocentric-transfer phase. 




Time, t, days 

(t>) Sensitivity of final angular velocity to thrust- 
vector variations. 

Figure 9. - Continued. Thrust-vector weighting functions 
for heliocentric -transfer phase. 






Time, t, days 

(c) Sensitivity of final radial position to thrust- 
vector variations. 

Figure 9. - Continued. Thrust-vector weighting functions 
for heliocentric-transfer phase. 
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(d) Sensitivity of final angular position to thrust- 
vector variations. 

Figure 9. - Concluded. Thrust-vector weighting functions 
for heliocentric-transfer phase. 
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(b) Radial-velocity time history. 

Figure 10. - Continued. Characteristics of capture 
spiral reference trajectory. 





Time, t, days 

(d) Radial-position time history. 

Figure 10. - Continued. Characteristics of capture-spiral reference 
trajectory. 
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(e) Angular -posit ion time history. 

Figure 10. - Concluded. Characteristics of capture -spiral 
reference trajectory. 
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(c) Sensitivity of final angular velocity to thrust- 
magnitude variations. 

Figure 11. - Continued. Thrust-vector weighting functions 
for capture-spiral phase. 
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Figure 12. - Continued. Simulated escape guidance maneuver^ thrust- 
magnitude perturbation, BF, 2.32xl0“ 3 newton. 





Final-radial-position error, 6r 
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(d) Final-radial-position characteristic. 


Figure 12. - Concluded. Simulated escape guidance maneuver; 
magnitude perturbation, SF, 2.32x10“^ newton. 
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